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ABSTRACT 



The aim of this work is to investigate properties of solutions of Fokker - Planck 
equation in the context of continuum mechanics. We show that average quantities, calcu- 
lated for these solutions approximately satisfy equations of isothermal motion of viscous 
ideal gas. 
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1. Introduction 

The statistical mechanics of systems described by kinetic equations and generalized Fokker - Planck 
equations is currently subject of active research (see [1]). 

In this work we study rather special kind of medium without interaction of particles. To describe the 
state of this medium we use density distribution function n(x h v,) in the space of Cartesian coordinates x t 
and corresponding velocities v, . 

We define average of any function / of these variables by 



Defined in such a way average value, depends obviously only on spatial coordinates x,. Macroscopic 
observer deals with these average values and traits them from the continuum mechanics field theory point 
of view. He perceives average / as volume density of some additive function F 
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The most interesting functions are functions, which satisfy some conservation laws. That is time de- 
rivative of F is a sum of volume and surface integrals 

~dt = \ \ \ d x idx2dxT, + J" J"^ Flidx 2 dx 3 + Fl 2 dx 3 dxi + Fl i dx l dx 2 - (3) 

D 

where Fl { is the current of /. 

Our aim is to construct the set of conservation laws for special case, when distribution function 
n(x,, v,) satisfies Fokker - Planck equation. 



2. General definitions of macroscopic parameters 

The density of mass is the average of unity 

1 = p. (4) 

The current of mass vector is 

v~i = pui. (5) 

where m, - average velocity. 

(5) represents at the same time the vector momentum density. 
The current of momentum tensor is 

J H = WJ- (6) 

Decomposition of velocity gives 

V; = W; + w,-; (7) 

where w i - chaotic component of velocity. Macroscopic observer considers chaotic motion as heat. The 
average of chaotic component w] is zero. 

Proceeding with (6), we have 



Jy = (Uj + Wj)(Uj + w j) = piijUj + WiWj. (8) 
The tensor of stresses <7y is defined by 

^ = -WJWJ; (9) 

- recall classic kinetic gas theory formula for pressure. 

Using this definition, we can write expression for the momentum current tensor in the final form 

= pUjUj - <7jj. (10) 

From macroscopical point of view the tensor of stresses describes the interaction between portions of 
continuum. It is remarkable, that though interaction between particles is missing, macroscopic cells are 
interacting nevertheless. This interaction is due to chaotic momentum transfer from one cell to another. For 
interacting particles additional term should arise in (9). 



The volume density of energy is 



e=^v k v k ; (11) 



because the only form of energy for non-interacting particles is their kinetic energy. 
Decomposition (7) of velocities leads to decomposition of full energy 



1 1 

e — — pu k u k + — w k w k = K + E. (12) 

The first term in (12) is kinetic energy density 

K=^pu k u k . (13) 

The second term in (12) is energy of chaotic movement or internal energy volume density 

E=^w k ~W k -. (14) 

For our simple model of non-interacting particles following relation between density of internal 
energy and pressure is valid 

1 2 

P = -j(7kk = jE; (15) 



where p - the hydrostatic pressure. 



The current of energy is 



1 1 1 1 

Fi = - v k v k Vi = - pu k u k Ui + - w k w k ui + u k w k Wi + - w k w k Wi = (16) 

1 

= (K + E) Uj - u k a ki + - w k w k Wj. 



3. Expressions for macroscopic parameters through Fourier transform of density 

Let us denote by M the Fourier transform of density 



oo oo oo 



M(t,x h qj)= — ^-3 J J | nexpH v k q k ) dv^dv^dv^ (17) 



-oo — oo — oo 



where q k - velocities momentum variables. 

Using this definition, we can express macroscopic parameters in terms of Fourier transform of den- 
sity and its derivatives: 

- density 



oo oo oo 



Af(f,x 1 ,x 2 ,x 3 ,0,0,0)=^3 \ \ \ nd Vl dv 2 dv 3 = ^ ; (18) 

— oo —oo —oo 



average velocity 



dM -i f f f -ipu k 

_(,,X 1 ,* 2 ,JC 3 ,0,0,0) = ^3 j J j Wt rf, lrf v 2 rfv3 = ; (19) 

— oo — oo -co 

current of momentum tensor 

_o oo oo oo 

d M -1 f f f 
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^ * 2 , * 3 , 0, 0, 0) = j j j ^,^3 = _i ; (20) 

volume density of energy 

9 oo oo oo 

d M — 1 f f f ~ e 

_(,,x I ,x 2> x 3 ,0,0,0)=^3 j j j mw^Mv 3 = ^ 



-co -co -co 



current of energy is 

-1 CO CO CO 

d 3 M i r r r 2iF,- 



— oo — oo — oo 



4. Fokker - Planck equation and conservation laws 

Fokker - Planck equation with damping force is: 



where a - damping coefficient. 

Multiplying (23) by exp(-z v k q k ) and integrating over velocities, we obtain 



(21) 



d(i d(hdqk (*,z 1 ,* 4 ,*3,0,0,0) = ^ j J" jn Vi v k v k d Vl dv 2 dv 3 =^-L. (22) 



3n dn dn d 2 n 

— + v t a - 3 an = - — - — ; (23) 

at ax k dv k ov m ov m 



dM 

~df + i 


( d 2 M d 2 M d 2 M > 


+ a 


?1 dq~ + qi dq~ + qi dq~)~ k \ qi+q2 + qi 


M. (24) 


ydx^dqi dx 2 dq 2 dx^dq^ j 


Substitute q k = to (24), use (18-19) and get 






dp | 3(pM t ) _ 
3f 3x t 


(25) 



This is the mass conservation law. 



Differentiate (24) by q h substitute q k = 0, use (19-20) and get 



d(pui) dJij 

+ — - + apu t = 0. (26) 



dt dxj 
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Use (10) and get 



d{pu t ) djpUjUj) day 

■H - — - — ^- J - + apu i -0. (27) 



dt dxj dxj 



This is momentum conservation equation. To obtain another form of equation substitute dp/dt from 



(12) to (13) and get 




diij 


+pu Jdx~ 


- -z-^- + apu; = 0. 


(28) 



Differentiate (24) by qj,q k , substitute q k =0 and get 

dM ik dM nk 

+ i + 2a M ik = -2k 5 jk M. (29) 

dt dXi 

where My... - derivative of M on q h q j, ... in zero point. 

Contract (29) on j, k and use (22-24) 

1 dJ kk dFj 

^ !" ^ a Jkk = _ 3fe p. (30) 

2 dt dx t 



de dF t 

— + ^-2ae = -3kp. (31) 
dt dXj 



Use expression (16) for F t 

d(K + E) d ( 1 ^ 

— dt — + dx~r K + E)u> ~ UkCTki + 2 WftWW ~ 2a ( K + ^ J = ~ 3k P ' 



(32) 



Contract equation (27) with n,- and get 



d(pu f ) d(pu iU j) day 

m, — Vui — - — - — Ui — - + apu t Ui = 0. (33) 

dt dXj dxj 



5. Equation of state 

The system of mass conservation law and momentum conservation equation is not closed. To close 
the system, we need express stresses through density, average velocities and their derivatives - equation of 
state. 



For this purpose differentiate (24) on q ( and q j and substitute q k = 0. 



dM,; dM uk 
dt dx k 



where My - M derivatives in zero point; Sy - Kronecker delta. 
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Equation (34) contains the third derivative. If local velocities distribution is normal (Maxwell with 
nonzero average), then Fourier transform of this distribution is also normal 

M = ex.p\^A ij q i qj + B k q k + c\ (35) 



We take (35) as working hypothesis. To give it some justification, we refer to fundamental solution 
[1] of Fokker - Planck equation. It is obvious from physical considerations, that after sufficiently long time 
of evolution solution for arbitrary initial conditions asymptotically converges to fundamental solution with 
the same centroid and momentum. 

If (35) holds, we get after simple calculations 

M ijk = \A ij B k + A jk B i + A ki Bj + B^jB Jexp(C). (36) 



Substitute (36) to (34) and get 



3<t j; 3<r,: du k du: dui 

~dF + Uk dx~ k +(Tij dx~ k + "* M + *> dx~ k + 2aa « = - 2kpS '- (37) 



To obtain (37), we used (25), (27) and (28). 

To simplify (37), we consider the static solution. In this case all velocities and their derivatives are 
zero and stresses are equal to hydrostatic pressure a t j = -{kla)pSij - see below section 'Static solution of 
Fokker - Planck equation'. We consider approximately static solutions and use linearization method to 
obtain state equation for this case. 

k 

<Tij=-- pSij + Scjy. (38) 
a 



Linearized form of (37) is 



l J-+2aSa lJ = p- — L + --U (39) 



dt 11 a \dxj dx/ 



We used the fact, that group of terms (k/a)[(dp)/(dt) + u k (dp)/(dx k ) + p(du k )/(dx k )] vanishes in con- 
sequence of the mass conservation law (25). 

Equation (39) is ordinary differential equation with constant coefficients and variable righthand side. 
If righthand side changes sufficiently slow, solution of differential equation exponentially converges to 

k 1 ( diij diij 



** = '*2 (40) 



Join (40) with (38) and find state equation 
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k k 1 ( du: du i , 



a 2 2 J 



We see, that (40) coincides with state equation of compressible viscous fluid. Relation between den- 
sity and hydrostatic pressure is p = (lc/a) p . 

It coincides with ideal gas isotherm. It was expected, because Fokker - Planck continuum is in ther- 
mal equilibrium with enclosing motionless medium. 

Kinematic viscosity is equal to 



k 






(42) 







6. Static solution of Fokker - Planck equation 

In this section we briefly discuss the simple and clear case of static solution. Namely we calculate 
average quantities for uniform static solution of Fokker Planck equation 



n = n Q e 2k - ™. (43) 
This is Maxwell distribution with zero average (see (47) below). Its Fourier transform is 

k 

M = M (t e~2a q "' q " '; (44) 



where 



x3/2 

M = n \—\ . (45) 



Density is equal to 



2.T)\W (1 »„ ^ J . (46) 



Average velocities are zero 



«* = 0. (47) 



The current of momentum tensor is 

k k 



Jy = (2xfM -S,j = p- S,j- (48) 



a a 
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and as (49) holds 



pu t uj = 0; 



(49) 



stresses are equal to 



-(2;r) 3 M - Sy =-p- Sy. (50) 
a a 



Hydrostatic pressure is 

k k 



p = (2x) 3 M Q - = p - . (51) 
a a 



(51) is the equation of state and we easily identify it with equation of state of ideal gas. Comparison 
of (51) with original equation of state of ideal gas gives known relation 

- = k B T; (52) 
a 

where 

k B - Boltzmann's constant; 
T - absolute temperature. 



The volume density of energy is 



k k 

e = 3/2(2^) 3 M - = 3/2p - = 3/2pk B T; (53) 
a a 

The kinetic energy density is equal to zero, because average velocities are zero. Therefore internal 
energy volume density is 

k k 

E = 3/2(2^) 3 M - = 3/2p - = 3/2pk B T; (54) 
a a 
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7. Fundamental solution of Fokker - Planck equation 

In this section we discuss macroscopic properties of fundamental solution (see [2]). Fourier transform 
of fundamental solution on both space coordinates and velocities is 



N = — -2 exp <^ —i 



I Pi , -at/ Pi \ 

XoiPi + Voi — +e (qi ) 

a a 



x exp < -k 



- 2 p iPi + 4(1" e~ at ) Pi ( qi - ^) + ^ (1 - e- 2 "<) ( 9j - (q, - 
a 1 a 1 a 2a a a 



(55) 



where p t - variables, associated with space coordinates; 

q t - variables, associated with velocities; 

and 



e = 2at-(l-e-°") (3-e" af )- 



(56) 



To obtain expression for Fourier transform M only on velocities, we need another set of variables 



x i = Xi -( XQi + V ^(l-e- at )). 
a 



(57) 



Fourier transform M of fundamental solution on velocities is 



M ■■ 



1 / 3 V>' 

1 ( a \ 



W\&ke) expHe " V0iqi)X 



(58) 



xexp 



\ (it ViV * 1 ia{ 1 '" " v + \ I a/ ( 1 - c 2a ' » - 2 n - ,: yn<i> 



Derivatives of M in zero point are 



1 / 3 V- 

1(a) 



(27rf \2nk0 ) CXP { 2k 



1 a" 



(59) 



■ / 3 x 3/2 / i 3 \ 



M <=(2^2^J exp^--x rt j|v 0ie - + «(l- e -) 2 ^|. 



-1 f a 3 f /2 f I a 3 \ 

M « = w{^) exp {--o2k^r (61) 

f v 0i e- a > + a(l- e-'f | Yv 0;e - + «(1 - e""') 2 y W ^ («' d " "2(1- ^) 2 ] . 
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We get from (19) expression for density 

( 3 ~\ 12 ( 1 3 ^ 

p = {2n?M «\^ke) ^{-oYk^} (62) 

After sufficiently long time we can drop exponentially small terms in (62) and get approximate solu- 
tion 

6' = 2at. (56') 



x\ = x ; - (x 0i + — ). (57') 
a 



3 \ 3/2 / j 3 



We see, that (62') is fundamental solution of diffusion equation (when v , is independent from x ,). 
So for sufficiently long times density approximately satisfies diffusion equation. 



Average velocity vector is 



{ MX 1 , = v m e- a < + a(l-e- a 'f^. (63) 

-ip j e 



The first term in (63) is uniform translation. It's velocity exponentially tends to zero. The second term 
is dilatation. After sufficiently long time the center of dilatation is (57') and dilatation rate is 0(t) (see 
(56')). 

Stresses are 

.3 A 3/2 / , „3 \ 



o„ = pu,u , + (2x? M, = [ - ^ j [at{\ - e- 2 "') - 2 (1 - e- a 'f j ) exp(- * ^ x, .v ( | o , . (64) 



Strain rate tensor is 



^- 2 ^ + ^:J- e (1 - e (65) 

We see, that (41) holds, when we drop exponentially small terms. 

In order to check momentum conservation law, we perform following calculations 

dp_ + d(puj) = 
dt dx t 
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3 \ de \ dd a _ _ a _ dx k \ ( a _ _ a 2 x k 2 1 



where 



di 2a(l-e- a ') 2 ; (67) 



and 



dt 

This proves (27) for our special case. 



M ' --v 0i e- a >. (68) 



d(pui) d(pUiUj) day 

■H — - — — J - + apu i = (69) 



dt dxj dxj 

= p\ -ui iH e- at + 2a 2 (l - e-°") ^ e~ at - a{\ - e~ al f ^ e~ at - a(l - e~ at f ^ 2«(1 - e""') 2 |+ 
+p(v 0i e-<" +a(l- e-f |]a(l - e"*) 2 * - p(- ^ *,)(" ^)(«f(l - e" 2 "') - 2 (1 - e~ at f ] + 

+ap| v 0ie - a '+«(l- e - a ') 2 f 
^ fy (1 - e""') 4 -2y faf(l - e" 201 ') - 2 (1 - e""') 2 1+ a 2 (l - e~ a, f + 2a 2 (l - e"" -2^(1- e""') 4 



P f (- y (! - e " a ') 4 " 2 y I «'< 1 " ' ) - 2 ( 1 - c ) 2 i + a 2 < 1 - c- 2 " 1 ) I = 



2 

X; a 



p — — | -(1 - £>""') - 2at(l - e~ za ') +4(1- e""') + (1 - e~ za, )(2at - (1 - e^")(3 - e~ a! )) 





PJT {1 ~ e "" )2 l " (1 " e ~ a,)1 + + e ~ a '* 3 ~ e ~ a "> 1= °- 



2 



Another conservation laws we can check in the similar way. 
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DISCUSSION 

We obtain equations of motion of Fokker - Planck continuum. Local average velocities and stresses, calcu- 
lated for Fokker - Planck equation solutions, satisfy these equations. 

With use of non strict arguments we deduct, that state equation of Fokker - Planck continuum is approxi- 
mately state equations of isothermal viscous compressible fluid. 

We obtain relation between kinematic viscosity of this fluid and coefficients of Fokker - Planck equation - 
coefficient of damping and coefficient of diffusion in velocities space. 

In the last sections macroscopic parameters of Fokker-Planck flows for static solution and fundamental 
solution are calculated. 
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